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We consider the class of risk measures associated with optimized certainty equivalents. 
This class includes several popular examples, such as CV@R and the entropic risk measure. 
We develop numerical schemes for the computation of such risk measures using Fourier 
transform methods. This leads to a very competitive method for the calculation of CV@R in 
particular, which is comparable in computational time to the calculation of V@R. 

Introduction 

The quantification of risk is more than ever a central issue in modern asset management. The increasing 
volume and complexity of financial instruments has raised the need for coherent but also accurate and 
efficient risk measurement methods. In the banking industry, a vast amount of positions and portfolios 
have to be assessed daily, which makes the computational speed of risk measurement methods a matter 
of paramount importance. Starting with Value at Risk (V@R), the goal of risk measures was to quantify 
the minimal amount of capital required in order to recover from unexpected large losses. V@R became 
very popular - see also the Basel II capital requirements - and is nowadays a standard instrument in 
the industry mainly for two reasons. First, it has an apparently obvious financial interpretation: it is the 
minimum amount of capital to be added to a position so as to push the probability of losses below a given 
level. Second, it has an easy and fast implementation: given a portfolio distribution, it simply accounts 
for the computation of the quantile of this distribution at a given level. However, V@R has a very serious 
deficiency, namely, it does not fulfill the basic property of diversification one should expect from a risk 
quantification instrument. Indeed, it may well happen that V@R delivers a lower risk for a portfolio 
concentrated in a single asset rather than one diversified into several assets. 

In order to overcome this drawback, Artzner et al. [1] introduced an axiomatic approach to coherent 
risk measures inciting diversification. An important example of such a risk measure is the Conditional 
Value at Risk (CV@R), which is strongly related to the Average Value at Risk (AV@R) and the Expected 
Shortfall. The seminal paper on coherent risk measures [1] was later generalised to monetary convex risk 
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measures by Follmer and Schied [16] and Frittelli and Rosazza Gianin [18], providing new examples, 
the most prominent of which is the entropic risk measure. This new class of risk measures attracted 
increasing attention from the financial industry as examples and numerical methods were studied, see, for 
example, Rockafellar and Uryasev [20], Dunkel and Weber [13] and Bardou et al. [2]. 

Compared to V@R though, coherent and convex risk measures are often more difficult to implement 
and costly in terms of computational time. Taking the example of CV@R, rather than looking for the 
quantile of the distribution at one level, it accounts for an integration of the quantile function over an 
interval, increasing thereby the numerical complexity. 

The goal of this paper is to focus on a class of risk measures - the optimized certainty equivalents - 
introduced by Ben-Tal and Teboulle [4, 3], and use Fourier methods to compute them efficiently. The 
first reason for the choice of this class is that it entails most of the classical examples - CV@R, entropic 
risk measure, and monotone Mean- Variance among others. The second reason is that due to its nice 
smoothness properties it provides a fairly easy scheme for computation which can be resumed in the 
following two steps: 

• Solve an allocation problem by use of a one dimensional root finding method; 

• Based on this optimal allocation, compute an expectation. 

Given these two steps and the explicit analytical formulas related to them, we make use of Fourier 
transform methods in order to compute these risk measures efficiently. Transform methods have been in- 
troduced to mathematical finance for option pricing, see Carr and Madan [5] and Raible [19] for instance, 
and have proved a very efficient tool when the Fourier transform of the underlying probability distribution 
is known. This is the case, in particular, for infinitely divisible distributions (i.e. Levy models) and affine 
processes. In the context of risk measurement, the application of transform methods has been largely un- 
explored. For this class of risk measures, they turn out to be a very efficient tool for computations as well. 
In particular, the calculation of CV@R using Fourier methods has similar computational complexity as 
the computation of V@R, thus both risk measures can be computed in almost the same amount of time. 
This should be a further argument supporting the use of CV@R in practical applications. 

This paper is organized as follows: in section 1 we discuss optimized certainty equivalents and their 
connections to risk measures. In section 2, we discuss the computation of optimized certainty equivalents 
using Fourier methods, and deterministic root-finding algorithms. We concentrate on the entropic risk 
measure, conditional value at risk and polynomial risk measures. 

1 Optimized Certainty Equivalent 

Let (CI, JP, P) be a probability space. By L° we denote the set of random variables identified when they 
coincide P-almost surely. By L p we denote the set of those random variables in L° with finite p-norm. 

Definition 1.1. A function I : M — > R is called a loss function if 

(i) I is increasing and convex; 

(ii) 1(0) = Oandl e dl(0); 

(iii) l(x) >bx + c and l(x) > b'x + d for all x e K for b > 1 > b' and c, c' < 0. 
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By I* we denote the convex conjugate of I, that is, l*(y) = s\rp xeR {xy — l(x)}. 

Let ip : M + — > K + be a convex function such that ip(0) = and f(x) > l(x) for all x G R + . In 
particular, <^ is a Young function and we denote by <p*{y) = swp x>0 {xy — (fi(x)} its convex conjugate. 
Following [9, 8], we define the Orlicz heart 

X v := {X G L° : E [ip{c \X\)] < +oo for all c> 0} (1.1) 

which is for the P-almost sure ordering and the (^-Luxembourg norm 

<lj (1-2) 

a Banach lattice. The norm dual of X v is the Orlicz space 

X* := {Y G L° :E[ip*(c\Y\)] < +00 for some c> 0} (1.3) 

with the Orlicz norm 

\\Y\\;:= S up{E[YX]:\\X\\ v <l}, (1.4) 

which is equivalent to the Luxembourg norm ||-|| .. Since (f(x) > l(x) > x for all x 6 M + , it follows 
that X v C L 1 . 

We denote with (P) the set of those probability measures on which are absolutely continuous 

wit h respect to P and whose densities are in X*. We consider risk measures in the following sense. 

Definition 1.2. A risk measure is a function p : X v — > [—00, +00] which is 

(i) quasi-convex: p(XX + (1 - X)Y) < max{p(X), p(Y)} for all X, Y e X v and A e]0, 1[; 

(ii) monotone: p(X) > p(Y), whenever X < Y for X, Y G X v . 
A risk measure is called monetary if it is 

(iii) cash additive: p(X + m) = p(X) — m for all m G R and all Ig^. 

It is known that any monetary risk measure is automatically convex, see [11, 18, 7, 12] and the references 
therein. Monetary risk measures were introduced by Artzner et al. [1] and extended to the convex case by 
Follmer and Schied [16] and Frittelli and Rosazza Gianin [18]. 

Given a loss function I, we define the Optimized Certainty Equivalent (OCE) introduced by Ben-Tal 
and Teboulle [3, 4] - to which we refer for further interpretation - as follows 

p(X):=ml{E[l(r ) -X)}-r ] }=m{ S l (r ] ,X), X e X v . (1.5) 

whereby 

Si(v,X):=E[l(r)-X)]-r], r, G R and X G X v . (1.6) 

The following proposition is known, up to minor differences in the assumptions. See [4, 6] for the 
case X v = L°°, [10] for the case where I is differentiable, and [9] for the computation of the dual 
representation in the general case. For the sake of readability, we provide a short proof based on results 
in [8]. 



\X\ 



inf { a > : E 



\X\ 
a 
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Proposition 1.3. Given a loss function I, the Optimized Certainty Equivalent is a lower semicontinuous 
cash additive risk measure taking values in R. Furthermore, for any X £ X v , there exists an optimal 
allocation rj* :~ i]*(X) £ R such that 

p(X):=E[l(ri*-X))-r l * (1.7) 

and this optimal allocation rj* belongs to [ess inf X, ess sup X] and satisfies 

E ( V * - X)] < 1 < E (r,* - X)] (1.8) 

where l'_ and l' + denote the left- and right-hand derivatives of I, respectively. 
Finally, it has the representation 

XeX. (1.9) 

If I is differentiable, this supremum is attained for Q* £ M-i.ip{P) where the density is given by 
dQ*/dP = V (77* - X). 

Proof. Since l(x) > x and X v C L 1 , it holds Si(t],X) > E[-X] > -oo, hence p(X) > -oo. On the 
other hand, Si(0,X) < E[l{X~)] < E[<p(\X\)] < +oo since X £ X v . Hence p(X) < +oo. 

Let us show that we have an optimal allocation determined by means of relation (1.8). Given X £ X v , 
the function r\ h-> Si(r), X) is real valued and convex. Furthermore it holds 

Si {rj, X) > E [-bX + br] + c]-r]>(b- l)rj - bE[X] + c 



P{X) 



QeM, v ,{P) 



max <^ E Q [-X] - E P 



dQ 
dP 



which goes to +oo as rj goes to +oo since b — 1 > 0. A similar argumentation with b', c' implies that 
Si(ri, X) goes to +oo as i] goes to — oo since b' — 1 < 0. Hence, there exists a minimum rj* £ R such 
that (1.7) holds. A straightforward argumentation shows that rf £ [cssinf X, esssupX]. Finally, this 
optimal allocation fulfills the first order optimality criteria 

Um S l ( V *+e,X)-S l ( V ,X) <Q<] . m S l (r ] *+e,X)-S l (r,,X) 

A straightforward application of Lebesgue's dominated convergence theorem allows to interchange limits 
and expectations and get relation (1.8). 

The fact that p is a cash additive risk measure is well-known, see [4]. The conditions of [8, Theorem 
2.2] are fulfilled and it holds 



P(X) 
a(Q) := 



max {E Q [-X] - a(Q)} , X £ X v 



(1.10) 



where 

sup {E Q [~X]-p(X)}, Q£M l , v .{P). 

However, since X v is a decomposable space in the sense of [21, Definition 14.59] and I is a normal 
integrand, we can apply [21, Thereom 14.60] yielding 



a{Q) — sup 





dQ Y 




7dP . 



V -E[l(-X + V )} 



sup ^ E 



dQ 



dP 



sup <^E 



rj>= sup < E 



dQ 
dP 



V [i-E 



dQ 
dP 



= E 
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which shows equation (1.9). As for the representation in terms of the optimal density, we refer to [10]. □ 



2 Numerical Computation of Optimized Certainty Equivalents 

The aim in this section is to develop numerical methods for the computation of optimized certainty equiv- 
alents based on transform methods and deterministic root finding algorithms. In general, the computation 
of the optimal allocation -q* and the risk measure p(X) in Proposition 1.3 can be performed in two steps: 

Step 1: use a deterministic root finding algorithm to compute rf in (1.8), combined with transform 
methods for the computation of the expectations; 

Step 2: use transform methods once more to compute the expectation E[l(rf — X)} and thus p{X). 

The terminus "transform method" indicates any method that uses the characteristic or moment generating 
function of the random variable X for the computation of expectations of functions of X. This includes 
the Fourier transform method of Carr and Madan [5], the Laplace transform method of Raible [19] and 
the cosine series expansion of Fang and Oosterlee [15]. Similarly, the term "root finding algorithm" refers 
to any method for determining the root of a function; e.g. bisection, secant or Newton's method, cf. Stoer 
and Bulirsch [23] for an overview. 

We start by describing the first main tool used in this section, which is the Fourier representation of 
expected values of functions of random variables. This result, which generalizes the classical Parseval 
formula, dates back to the work on transform methods for option pricing and is particularly suitable 
for models with known characteristic function, e.g. Levy and affine models. We will follow the work of 
Eberlein et al. [14] closely, while we refer to Schmelzle [22] for a comprehensive overview and numerous 
references. 

Let I be a loss function as described in the previous section and denote by Ir the dampened loss 
function, defined by Ir(x) :~ e~ Rx l(x), for R S R. Moreover, let / denote the Fourier transform of a 
function /, i.e. f(u) — J e lux f(x)dx, and Mx the (extended) moment generating function of X, i.e. 
Mx(u) = E[e uX ], for suitable u G C. By L 1 , resp. L^, we denote the set of measurable functions on 
the real line which are integrable, resp. bounded, continuous and integrable, with respect to the Lebesgue 
measure. We also denote by 5° the interior of a set S and by the imaginary part of the complex 
number z. 

Theorem 2.1. Let X e L° and define 

1 := {R G M : M X (R) < oo} (2.1) 
J := ji? € R : Ir G L\ c andT R E L x \ . (2.2) 

Under the assumption 

(A-I) IflJ/fl, 
the following representation holds true: 

E[l(r) - X)] = i J 4 R - lu ^M x {iu - R)T{u + iR)du. (2.3) 

Proof. See [14, Thm. 2.2]. □ 

Remark 2.2. Consider a random variable X such that G 1° and a loss function / that has polynomial 
growth. Then, it follows directly that X G X v . ♦ 
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2.1 Entropic Risk Measure 

We start by considering the entropic loss function 

e fX - 1 

with 7 > 0. The derivative and the conjugate function are 

y My) y - 1 



l'(x)=e< x and Z*(j/) = 



7 7 



In this case, the optimal allocation rf and the risk measure p{X) can be computed explicitly and are 
provided by 

V* = - 1 -ln(E[e-^]), 
p{X) = -r ] * = Un(E [e-^]) 

'dQ 



sup <; £ Q l-X] - E Q 
QeXi,„»(P) 



In 



Note that in Levy and affine models the quantity E[e ~ (X ] is known explicitly in terms of the Levy- 
Khintchine formula or the affine transform formula, hence also rj* and p(X) are computed explicitly. 

2.2 Conditional Value at Risk 

The most interesting example from the point of view of practical applications is Conditional Value at 
Risk, also known as Average Value at Risk or Expected Shortfall. These notions coincide if X has a 
continuous distribution, see Follmer and Schied [17, Corollary 4.49]. The loss function in this case is 

, „ I 71X if x < 

l(x) = - 7l x" + 7 2 z + = <^ " (2.4) 

I 722; if x > 0, 

with 72 > 1 > 71 > 0. The left-hand derivative equals 

l'-(x) = 7i 1 {z<o} +J2 1 {x>0}, (2-5) 

while the conjugate function is 

,»/ \ f° if 72 < x < 71 

t (s) = < (2.6) 
l+oo otherwise. 

In case 71 = 0, this corresponds to the standard CV@R with parameter I/72; see e.g. Rockafellar and 
Uryasev [20]. 

The optimal allocation can be computed explicitly, in terms of the quantile function q x of X, and is 
provided by 

V* = ix . (2.7) 

\72 - 71/ 
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The following representation for this risk measure is also standard in the literature 

1—71 

11— 71 1 

p(X) = --f 2 J qx{s)ds-j 1 J q x (s)ds (2.8) 

l-Tl 
11 -71 

= sup {Eq [-X] : 7i < dQ/dP < l2 ] . (2.9) 

QeM ltV *(P) 

In particular, for the special case of CV@R with parameter A = I/72, it holds that rf = qx (A) and 

A 

CV@R X {X) = -\ [ q x (s)ds = sup \E Q [-X] : dQ/dP < M (2.10) 

The aim of the next result is to provide an alternative representation for p{X), using transform methods. 

Proposition 2.3. Assume that the optimal allocation 77* is computed by (2.7). Let X G L 1 be a random 
variable such that I 7^ with G 1°. Then, p admits the following representation 

= t~ / r^^2 M ^( m - R ^ du - t~ / ; ^-nv M x( iu - R2)du - v*, (2.ii) 

2tt J {u + iRiY 2-k J {U + 1R2Y 

R R 

where R\ G Xf] (—00, 0) a«ii ]?2 6lfl (0, +00). /« particular, for CV@R we get 
CV^(X) = —j^ TW M xi iu-R) d u 

R 

w/iere A = I/72 R & T D (—00, 0). 

Proof. Since 77* is already computed using the first order condition (1.8) for the loss function (2.4), see 
(2.7), we will apply Theorem 2.1 directly to representation (1.7). We get 

p(X) = E [l( v < - X)] -r,*=E [- 7l (?f - X)- + 72(77* - X)+] - rf 
= - 7 i£ [(X - 77*)+] + l2 E [(77* - X)+] - rf 

= -7T I e (Rl ~ iu)r, *M x (iu- Ri)T 1 (u + iR 1 )du 
2n J 

7. 

+ £ / e(fe " m) "* M ^(« M - R2%(u + iR 2 )du - 77*, (2.12) 

R 

for i?i Gin Ji, i?2 Gin where we define the functions 

/i(a;) = (~x) + and Z 2 (a:) = + . (2.13) 

Now, we just have to compute the Fourier transforms of the functions l\ and Z 2 — which are reminiscent 
of the (payoff) functions corresponding to a put and a call option — , determine the sets J\ and J2 and 
show that the prerequisites of Theorem 2.1 are satisfied. 

We have that, for z G C with G (— 00, 0), the Fourier transform of l\ is provided by 



u 

h{z) = J e izx {-x) + dx = - J c lzx xdx = ~y/ ZX 

R -00 



(iz) 



e izx 



-\, (2-14) 
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while for I2 we get the same formula, that is 

h{z) = -\, (2.15) 

where now z e C with ^s(z) e (0, +00). 

Let us turn our attention to the corresponding dampened payoff functions, i.e. 

l hRl (x)=e- R ^(-x)+ and l 2 , R2 (x) = e~ R * x (x)+ . (2.16) 

Clearly, for R\ < and R 2 > 0, these functions are bounded and continuous, while from (2. 14) it follows 
directly that h, Rl 1 h,R 2 E L 1 . A direct computation shows that also h >Rl , fa,R 2 6 Indeed, 



|2 _ /"l; . _ /„M2 , _ f „-2R lX Jl • 



||ii,Hill£»= / IMi = / e" 11 A = 4^3 < °°< < 2 - 17 ) 

R -00 1 

while the computation for I2.R is completely analogous. We can also examine the weak derivatives of 
h,R, h,R', we have that 

( e -Ri*(R lX -l) for a; > 
^i,«x(ar (2-18) 
[0 fora;<0 

and we can directly deduce that dl\.Ri £ £ 2 , while the same is true for dl2,R 2 - Thus, we get that 
h,R-L 1 h,R 2 belong to the Sobolev space 

H^R) = {g EL 2 :dg exists and dg G L 2 } (2.19) 

and using [14, Lem. 2.5] we can conclude that lx,R i: h,R 2 are integrable. Therefore, J\ = (—00, 0) and 
J 2 = (0,+oo). 

Finally, since 6 1°, we have that X n J\ 7^ and in J2 7^ 0, hence assumption (A-I) is satisfied. 
The result now follows by substituting (2.14)-(2.15) into (2.12). □ 

Remark 2.4 (Numerical Efficiency). The Fourier representation (2.1 1) for the Conditional Value at Risk 
is numerically more efficient than the standard one in (2.8). Indeed, the latter requires to solve an op- 
timization problem, i.e. the computation of the quantile qx, for every grid point used in the numerical 
quadrature, while the former requires to solve only one optimization problem for the computation of 77*. 
Assume that the grid for the numerical integration has size N, the computational effort for the solution 
of the optimization problem is Mq, while the computational effort for the numerical integration is Mj, 
where typically Mj Mo- Then, the total computation effort for the two methods compares as follows: 

TCE (Fourier) = Mo + Mi vs TCE (standard) = N ■ Mo + Mi. (2.20) 

This also reveals that the computation of CV@R is not significantly more time consuming that the com- 
putation of V@R, since the bulk of the computation is the optimization problem (i.e. quantile or V@R) 
and not the numerical integration. This should be an argument in favour of using CV@R in practical 
applications. ♦ 
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2.3 Polynomial Loss Function 

Next, we consider the polynomial loss function defined by 

([l + x] + Y-l 
l(x) = - - (2.21) 

7 

for 7 6 N, 7 > 1. The derivative equals 

I'(aO=([l + a;] + ) 7 ~\ (2.22) 



and the conjugate function is provided by 

(i - ^y- 1 - 72/ - i 



for y > 



i*(y) = { 1 ' (2.23) 

oo otherwise 

In this case, neither the optimal allocation nor the OCE can be computed explicitly, and one has to resort to 
numerical methods for both. In particular, we will combine Fourier transform methods with deterministic 
root-finding algorithms. 

Proposition 2.5. Let X G L 1 be a random variable such that G 1°. The optimal allocation is deter- 
mined numerically by the secant method, that is % — > rj* as k — > oo, where 

Vk+i = Vk - /(%) • -77—^ — r (2.24) 

JKVk) - f{Vk-l) 

with 

(-V — IV f e (iu-R)(l+, lk ) 

/(*) = ^ J Mx{iu - fi ) (iu-Rp du - 1 (Z25) 

R 

Once rf is determined, the polynomial loss function risk measure admits the following representation 

far-lV f e (iu-R)(l+r]*) 1 

PW = J M *(™ R ) {iu _ R) , + i du - - - rT, (2-26) 

R 

where Rein (-oo, 0). 

Proof. We start by computing the Fourier transform of the following function: 

£( x ) = [( a - x) + ] n (2.27) 
for a G R, n G N. Integrating by parts iteratively we get, for z G C with G (— oo, 0), that 



£( x ) = j c lzx p{x)dx = J c lzx {a - x) n dx 

R -oo 

a 

= (a-x) n +- / e izx (a-x) n - l dx 

IZ -oo IZ J 



=0 



77 1 77 



e lzx dx = ,. - .. e"". (2.28) 



(iz) n J (iz) n+1 
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Following the same argumentation as in the proof of Proposition 2.3, we can show that the dampened 
function Ir belongs to L^ c and has an integrable Fourier transform for all R 6 (—00, 0). In other words, 
J = (-oo,0). 

Now, consider the function 



- 1. (2.29) 



According to (1.8), the root of this function determines the optimal allocation corresponding to the poly- 
nomial loss function (2.21). This root can be determined by a standard root-finding algorithm, e.g. the 
secant method, which states that r)k — > rf as k — > 00, where 

Vk+i = Vk - /(%) • -jj-A — — r- (2.30) 

fiVk) - f(Vk-i) 

Applying Theorem 2. 1 to (2.29), using (2.28) with a = 1 + 7] and n = 7— 1, and recalling that ID J ^ 
since G 1° , yields the representation (2.25). 

Once the optimal allocation has been determined numerically, we just have to combine (1.7), (2.21), 
Theorem 2.1, and (2.29) with a = 1 + 77* and n = 7, and direct computations yield representation (2.26) 
for the risk measure corresponding to the polynomial loss function. □ 
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